Metallic ferromagnetism in the Hubbard model: A rigorous example 
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We present the first rigorous example of the Hubbard model in any dimensions which exhibits 
metallic ferromagnetism. The model is a genuine Hubbard model with short-range hopping and on- 
site Coulomb repulsion, and has multi single-electron bands. In the limit where the band gap and the 
Coulomb repulsion become infinite, we prove that the ground states are completely ferromagnetic 
and at the same time conducting. 
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It has been believed since Heisenberg [l| that ferro- 
magnetism observed in nature is generated by quantum 
effects and Coulomb interaction between electrons. It is 
a challenging problem to confirm this scenario by show- 
ing that only short-range hopping of electrons and spin- 
independent Coulomb interaction can lead to ferromag- 
netism in the concrete setting of the Hubbard model [2| • 

Now many rigorous examples of ferromagnetism (or 
ferrimagnetism 3]) in the Hubbard model are known, 
and it is clear that certain versions of the model do 
generate ferromagnetism. An important class of exam- 
ples, now called flat-band ferromagnetism, was discov- 
ered by Mielke Q and then by Tasaki [B| . In these mod- 
els electrons occupy the lowest dispersionless band, and 
infinitcsimally small Coulomb interaction can lead to a 
complete ferromagnetism. Related models were found in 
Although the flat-band models are singular in the 
sense that the single-electron ground states have huge 
degeneracy, the mechanism which generates ferromag- 
netism is believed to be robust and physically realistic. 
Indeed the existence of ferromagnetism has been proved 
rigorously in related nonsingular models 0, S] • 

A common feature of all these rigorous examples of fer- 
romagnetism is that they describe insulators [9( . Metallic 
ferromagnetism, in which same electrons contribute both 
to magnetism and conduction, is clearly more interesting 
and challenging. As far as we know the only rigorous ex- 
ample of metallic ferromagnetism in the Hubbard model 
is that by Tanaka and Idogaki (lp| . who treated a quasi 
one-dimensional model using the Perron-Frobenius argu- 
ment 0. But the physics of one-dimensional electron 
systems is very special, and it is highly desirable to have 
examples in higher dimensions. 

In this Letter, we present the first rigorous example 
of metallic ferromagnetism in a version of the Hubbard 
model in any dimensions. The mechanism of ferromag- 
netism in the present model is basically the same as that 
in the previous models, namely, when one represents the 
system using a moderately localized basis, the Coulomb 
repulsion (in real space) generates both a repulsive inter- 
action and a ferromagnetic exchange interaction. Our 
model is a variant of the models in [!, [H, 0] and has 



multi single-electron bands, among which the lowest two 
mainly contribute to low energy physics (especially in the 
large band gap limit that we take). In the ground states 
the lowest band is half filled and exhibits ferromagnetism 
as in 0, H, 0| ■ The electrons in the second lowest band, 
which is partially filled, are movable and are coupled fcr- 
romagnetically to the electrons in the lowest band. This 
gives rise to ground states which are ferromagnetic and 
at the same time conducting. 

Although the basic mechanisms are similar, the math- 
ematical methods developed for the insulating models 
[ESQ never apply to conducting systems [H[ . We here 
develop a completely different variational argument. 

Definitions and main results: Let A be the d- 
dimensional L x • • • x L hypercubic lattice (where L is 
even) with unit lattice spacing and periodic boundary 
conditions. Let N = {(1/2, 0, . . . , 0), . . . , (0, . . . , 0,1/2)} 
be the set of d vectors of length 1/2 pointing in the 
positive direction of each axis. Then B = {x + 8 \ x € 
A, 5 E J\f} can be regarded as the set of mid-points of 
bonds in A. We construct a Hubbard model on the lat- 
tice A = A x {1,2, 3} U B x {1, 2} (where the triplicated 
lattice A x {1, 2, 3} consists of pairs (x, i) with x S A and 
i = 1, 2, 3, and the duplicated lattice B x {1, 2} consists 
of pairs (w, i) with w € B and i — 1,2). See Fig. [I] With 
each site z € A and spin index a =t, J,, we associate the 
standard fermion operator c Zi<T . 

It is convenient to define some fermion operators by 
combining the basic operator c ZjCr . For each x £ A, 5 € 
Af, i = 1, 2 and a =|, |, we define 
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where v > 0, is a model parameter, whose value does 
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FIG. 1: The lattice A for d = 1. Integer sites are triplicated, 
and half-odd-integer sites are duplicated. 



not play essential roles in the present work. See Fig. [5] 
These operators are designed in such a way that any 
electronic state on A can be written by a combination 
of a, b, and d operators. Moreover one has {a\b] — 
{a\d} — {b\d} = for any combinations of indices, 
and {a\ a , a VlT } = {b\ t(T , b y ^ T } = 5 X!V 5 aiT for any x, y G A 
and a,T =t,|- Note that, unlike in our previous mod- 
els 0, 0] , the a-operators satisfy the standard canonical 
anticommutation relations. 
We define Hamiltonian H by 



H = 



y ] ( sa x,o- a y,& t b x,a b y,cr) 
x,y(zA; \x — y\ — 1 
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with 



<r=T.J. 

= rt 



Note that ([5]) defines a gen- 



uine Hubbard model with short ranged (but admit- 
tedly complicated) hopping amplitudes. The model has 
several bands; the a-band with the dispersion relation 
e a (k) = — 2s X^i=i cos(fcj), the 6-band with e&(k) = 
v — It J2i=i cos(fcj), and the d-bands with higher ener- 
gies. We fix the total electron number to N c . 

Theorem — Let d = 1, 2, 3, ... be arbitrary and suppose 
that |A| < iV < 2lA| El and v > 2d(\s\ + 2\t\). In 
the limit u, U — > 00 [L3J, the ground states of ([5|) exhibit 
saturated ferromagnetism in the sense that they have the 
maximum possible total spin Stot = N e /2. 

One may replace the lower bound for v by better values 
which depend on the electron number. For example it is 
enough to have v > 2d (\s\+\t\) when 3|A|/2 < iV c < 2|A|. 

The electron number iV c = |A| corresponds to the half- 
filling of the lowest a-band, and iV c = 2|A| to the half- 
filling of the a and 6-bands. Therefore when t ^ and the 
electron number satisfies |A| < N c < 2|A|, the ferromag- 
netic ground states, which are indeed Slater determinant 
states, are conducting states with N c — |A| conducting 
electrons (or 2|A| — iV c holes) in the 6-band. 

Finite energy states: We shall describe a complete 
proof of the theorem. We say that $ is a finite energy 
state if (<£>,i/<I>) < 00 in the limit u,U — > 00. A finite 
energy state cannot contain any of the d-states since a 
d-electron costs an energy proportional to u, which be- 
comes infinite. Furthermore since U — > 00, a finite en- 
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FIG. 2: Components of the states corresponding to the spe- 
cial fermion operators indexed by x 6 A and w 6 B. We 
omitted the normalization factors for the a and b operators. 
Each state is localized within the unit lattice spacing. 



ergy state $ must satisfy for any z G A the condition 
'T-z,T n z,i c f > = and hence 

Cz, iCz , T $ = 0. (6) 

Let $0 be the state with no electrons. A compu- 
tation shows that C( x ^ ^C( x .3),f(- ■ -) a l | a l x^o = (3 + 
4aV 2 ) _1 (- • -)$0! where (• ■ •) is an arbitrary product of 
and 6^ except for a x *, a x , . This means that any 
state $ which contains a term with a x ^a x , cannot sat- 
isfy C( a . j3 ) xC( 2:i 3),| ( I ) = 0. Thus a finite energy state i> 
has no terms with a x ^a x ^. Likewise [lij ] we can show 
that 4> has no terms with b x ^b x , . Therefore any finite 
energy state is in the subspace 7i HC with the "hard core 
condition" , which is spanned by the basis states 

*(A *;B,t) = (J[ al Mx) ) ( J] bl <x) ) $ , (7) 



x£B 



where A, B are arbitrary subsets of A such that |^4| + 
\B\ = N c , and er = (a(x)) xeA , T = (t(x)) x( z B are arbi- 
trary spin configurations with a(x),r(x) G {|, 1} [l5| . 

For a state <!> to satisfy (O, it is not enough that 
$ G 7Y HC . By imposing © for other sites, we find that 
a finite energy state $ must satisfy the following local 
ferromagnetic conditions. When we expand $ as 



$= J2 *KA,a;B,T)*(A,a;B,T), 



(8) 



A,CT,B,T 



the coefficients tfj(A,cr: B,t) must satisfy 
ip(A,cr;B,T) — %jj{A,cr x ^y]B,T) for any 1,1/6 4 such 
that |a; — y\ = 1, and tp(A, tr; B, t) — ip(A, a x ; B, t x ) for 
any x G A n B. Here <r x ^ y is the configuration obtained 
from er by switching o~{x) and cr(y) in a '. Similarly 

cr x ,T x are obtained by switching cr(a;) and t(x) in 

cr,r [16|. These conditions are equivalent to infinitely 
large ferromagnetic couplings between neighboring a- 
electrons, and between the a-electron and the b-electron 
sharing a same site x. 

By 7i LF we denote the subspace of 7i HC consisting of 
states which satisfy the local ferromagnetic conditions. 
Note that for any <!> G 7Y LF the expectation value of H 
satisfies (<I>,7d$) = (<f>,i/ off $) with 
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Variational estimates: So far all of the arguments are 
straightforward variations of those developed for the sim- 
ple flat-band models Let us now turn to variational 
estimates, which are essential to our treatment of con- 
ducting states. 

Note that the above stated local ferromagnetic condi- 
tions relate the coefficients il)(A, cr; B, r) with common 
A and B. We can thus decompose 7i LF into a direct sum 
as H LF = ®^Na=N c -\A\^Na- Here i s the intersec- 
tion of 7i LF and the space spanned by the basis states 
*5>{A, a; B, r) with any A such that \A\ — N a , and arbi- 
trary cr, B, and r. Since the effective Hamiltonian (J9j) 
leaves the number of a-electrons invariant, we can deter- 
mine the ground state energy Eqq, variationally as 



East = 



min E(N a ) 

|A|<jV„<|A| 



E(N a )= inf ($,ff cff $) 

(10) 

When N a = |A|, a-electrons fill the entire A, and are 
coupled ferromagnetically. Since all o-electrons are cou- 
pled ferromagnetically to the a-electrons, we see that 
any state in Hj^j has the maximum possible total spin 
Stot — N c /2. It is also easy to see that J5(|A|) gives 
the lowest energy among the ferromagnetic states. In 
what follows, we shall prove that E(N a ) > -E(|A|) for 
any N a < |A|. This shows that the ground states have 
the maximum total spin, and proves our theorem. 

Let N a < |A|. We first note that on the space H.\f, 

H cS >H + (N c -N a )v- 2d\s\ (|A| - N a ), (11) 
where 



H = -t 



E 



u x,<j u y,<?- 



(12) 



x,y£A; \x — y\ — 1 

<r=T4 



To get the lower bound (fTTj) . we noted that each hole 
(i.e., a site in A not occupied by an a-electron) has a 
kinetic energy not less than — 2d\s\ 171 ]. 

Since H does not act on a-electrons, we have 

inf = min inf (13) 

*ew^, F ^ca <s,£hy 



ll*ll=i 



\ A \= N ° ||*||=i 



where 7i^ F is the intersection of TC LF and the space 
spanned by the basis states ^(A, cr; B, r) with the spec- 
ified A and arbitrary cr, B, and r. 

Note that, on Ti. 1 ^ (even on H HC ), we can bound H 

as 



H>H A -2d\t\(\A\-\A\), 
for any Ac A, where 



H A = -t £ 

x,y£A; \x—y\ = l 

°-=T4 



u x,cr u y,o 



(14) 



(15) 



is the hopping Hamiltonian restricted on A. To get the 
lower bound (fT4|) . we applied the bound | 2 CT (^x a^y," + 
h.c.) | < 1 (which is valid on 7i HC ) to all the hopping 
terms including any site in A\A. From fTOjl. (fTTjl . (fT3]) 
and (I14| . we have 

|A|)« + (|A|-JV a ){«-2d(M + |t|)} 
min inf ($,H A $). (16) 

\A\=N a m = 1 



E{N a ) > (N c 



We shall examine the infimum in (|16|) . Let us de- 
compose A into connected components as A = 1J™ =1 Ai. 
Within each Ai, all the a-electrons and o-electrons are 
coupled to have the maximum possible total spin be- 
cause of the local ferromagnetic conditions. Note that 
Ha allows 6-electrons to hop around only within each 
connected component Ai, and leaves those o-electrons on 
A\A unaffected. This means that the above ferromag- 
netic coupling is not disturbed by the application of H A . 
Therefore the infimum of the expectation value of H A 
taken over all states in 7i^ F can be estimated simply in 
the subspace spanned by up-spin electrons [ijj]. At this 
stage we can forget about the a-electrons, which have no 
kinetic energies in H A , and consider only the (now fully 
polarized) o-electrons. This leads us to 



inf ■ 

ll*ll=i 



$,H A &) = inf (*,Ha&), (17) 



<S>£H 



- -i 



11*11 = 1 



where 7i\ N is the space spanned by states of the form 
(ELes w hh an arbitrary B C A such that \B\ — 

N. An inspection shows that, in the space T~t\ N -\ A \> 
the number of movable electrons (which are on A, and 
hence acted on by H A ) varies from iV m i n — iV c — |A| to 
N max = min{\A\,N e -\A\}. 

Since unmovable electrons (which are on A\A) 
are not affected by H A at all, we see that 

((IlxeB &£, T )*o,#a(IL 6 b'&£,t)*°> = whenever \An 
B\ =^ \An B'\. Therefore we can evaluate the infimum in 
the right-hand side of (|17p in each subspace with a fixed 
number of movable electrons. Since unmovable electrons 
have no contributions to expectation values of H A , we 
have 



<s>en 



inf ($, H A $) 



mm 



A,JV C -|A 
11*11 = 1 



W min <W<JV max 

" " 11*11=1 



inf ($,#a$>, 



T 

A.N 



(18) 

where Tt A N is the space spanned by states of the form 

dlxG-B ^ |)^o w hh an arbitrary B C A such that \B\ — 
N. 

Let H\ be the hopping Hamiltonian (| 1 5[) with A = A. 
Since H\ and H A are equivalent when restricted to the 
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subspace Ti \ N , we see that 

inf (®,H A ®) = inf (<J>,iTA$)> inf ($,H A &), 



T 

A,N 



ll*ll=i 



T 

A . N 



ll*ll=i 



where the inequality follows from 1-L A 
then find, from dTTJ) , ([18]), and HU), that 



A,JV 
*ll=l 



.1. V C ^0\. V 



(19) 
We 



inf ($, > 



mm 

N 

^<N<N B 



inf ($,if A $). (20) 
*e«! 



l A,N 



11*11=1 



Let ei < £2 < • • • < £|A| be the eigenvalues of the 
hopping Hamiltonian H\ (which is (|15p with A — A) in 
ascending order. Since the energy spectrum has a plus- 
minus symmetry, we see that eg < if £ < |A|/2 and 
e^>0if£>|A|/2. The infimum in the right-hand side 
of (f20l is nothing but the ground state energy of a spinless 
free fermion, and is equal to X^=i e ?- -By minimizing this 
over N, we see that [19( 

Wmin 

(right-hand side of (20])) > (|A| - |A|) e + ^ e £ , (21) 

e=i 

where e is if N min > |A|/2, and is £Ar min if AT min < \A\/2. 

As for the lowest energy E(A) of the ferromagnetic 
states, one has 



E(\A\) = (N c -\A\)v + 



|A| 

J2 e ^ 



(22) 



since one simply fills all the a-states and the lowest 
— |A| states of the 6-band with up-spin electrons to 
get the lowest energy. Note that the total energy of the 
fully filled a-band is vanishing since there are no diagonal 
terms in the hopping Hamiltonian of the a-electrons. 
Combining (JT6J), {20]), ([21]), and (22]), we finally get 



E(N a )>E(\A\) + (\A\-\A\){v-2d(\s\ + \t\) + e}, (23) 

which implies the desired bound E(N a ) > -E(|A|) if N a < 
\A\ and v > 2d(\s\ + \t\)-e. Since e > -2d|t|, we get the 
condition for v in the theorem. The improved condition 
is obtained by recalling that e = if N m i n > |A|/2. 
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Let (■ • ■) be an arbitrary product of and 
6^ for sites other than x. Then we have 

C(*A)a<W).T(--0&s lT &ia*° = (1/2)0 -0*0, and 

Ei=i,2 c (^,i),i c (^,i),T('-') a i,<T 6 l,T & l,i <1>0 = (• • -)«£, CT *o- 

These show the claim. 

We order the sites in A in an arbitrary but a fixed man- 
ner. Products are arranged according to this order. 
Let x, y £ A with \x — y\ = 1 and set w = (x+y)/2. If (■ ■ ■) 
is an arbitrary product of and 6^ except for 



,t 



%,V a l,l> wenavec (»,2),i c (»,2),T(' ' 'Wx^al^o = ^ 2 (3+ 
4d^ 2 ) _1 (- • -){^T,o-5j.,t — <5j, CT <^,-r}<I > o. This, when applied 
to the expansion (8]), yields the first condition on ip. Like- 
wise, if (■ ■ ■) is an arbitrary product of and for sites 



other than x, we have C( xA)>l C( x ^)^(- ■ -)a x ^b x ^<& = 
{2(3+4di/ 2 )}~ 1/2 (- ■ •){5 tiCT (5 i , r -(5 i , CT <5 T>T }$o, which leads 
to the second condition on ip. 
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Let $| be any normalized state in Tt A F which consists 
only of up-spin electrons. Let S x = a x ^a x ,i + b x ±b x ,-[ 
and S~ = Yl x eA- be the spin lowering operators. De- 
fine S(n) = I\ x L\a( S *T x H^^T* with ". = 0,1 



and rii = 0,1,2..., where 
collection of these n's, and 
= (normlization) S(n)$-j 



n is a short hand for a 
define a normalized state 
. Then the space 



ojLF 

n A 



spanned by the states of the form $ n . Since Ha and S(n) 
commute, {& n > , Ha&h) equals ($ T , Ha$i) if n' = n and 
is zero otherwise. Thus the claim follows. 
Suppose that iV m i n > |A|/2. Then since te > for 
£ > Y m i n , we see that the desired minimum is e ?- 
When iVmin < |A|/2, we can bound the minium from be- 
low by J2f=i° e e + (Ymax — Y min ) ejv min > J2e=i" Q + 
(|A| — |A|)ejv min , where we used iV m i n = N c — |A| and 
Y max < Y c -\A\. 



